Semi-discrete composite solitons in arrays of quadratically nonlinear waveguides 
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Abstract 

We demonstrate that an array of discrete waveguides on a slab substrate, both featuring the \ 2 nonlinearity, supports 
stable solitons composed of discrete and continuous components. Two classes of fundamental composite solitons are identified: 
ones consisting of a discrete fundamental-frequency (FF) component in the waveguide array, coupled to a continuous second- 
harmonic (SH) component in the slab waveguide, and solitons with an inverted FF/SH structure. Twisted bound states of 
the fundamental solitons are found too. In contrast with usual systems, the intersite-centered fundamental solitons and bound 
states with the twisted continuous components are stable, in an almost entire domain of their existence. 



PACS numbers: 

Quadratically nonlinear (x^) media, continuous or 
discrete, provide favorable conditions for the creation 
of optical solitons. The wave-vector mismatch and x^ 
coefficient are efficient control parameters in this con- 
text, and the solitons display a variety of features due 
to their "multicolor" character. Accordingly, a great ef- 
fort has been invested in the study of solitons in con- 
tinuous EU3 ( as reviewed in 0,0) and (quasi-)discrete 
H S II II IE G3 E3 X (2) media. In both cases, the 
solitons can find applications to all-optical switching 
[Hill 111 and light-beam steering. [H 

Waveguide arrays, i.e., one-dimensional (ID) discrete 
systems, exhibit properties that are absent in continuous 
media, such as anomalous or managed diffraction. |ll| Ac- 
cordingly, discrete solitons are drastically different from 
their counterparts in continuous media, as was first pre- 
dicted in the context of the x^ 3 ) nonlinearity Here 
we propose semi-discrete composite solitons in x^ op- 
tical systems, that contain both discrete and continu- 
ous components, each carrying either the fundamental- 
frequency (FF) or second-harmonic (SH) wave. Wc 
demonstrate that stable semi-discrete solitons can be 
readily formed in a waveguide array coupled to a slab 
waveguide, both structures being made of a quadrati- 
cally nonlinear material. We study two most interesting 
species of semi-discrete solitons. Type- A ones consist of a 
discrete FF component in the waveguide array, coupled 
to a continuous SH component in the slab waveguide. 
Conversely, type-B solitons feature continuous FF and 
discrete SH components in the slab and discrete array, 
respectively. 

The proposed setting is displayed in Fig. ^ It in- 
cludes the periodic array of waveguides, with a spacing 
xq, mounted on top of (or buried into) the slab waveg- 
uide. Both the array and slab are made of a x^ material, 
such as LiNbC>3 or KTiOP04. A rigorous coupled-mode 
theory for such composite waveguides can be developed 
by a straightforward generalization of that available for a 
single discrete waveguide coupled to the slab ^?|]- Thus, 



we arrive at a system including a set of ordinary differen- 
tial equations for the discrete array, coupled to a partial 
differential equation for the slab. In case A, the coupled 
equations take the following normalized form: 



i^r + Q (<Pn-i + <t>n+x) + <ft*(C, n) = 0, (fa) 



id 2 * 

'Ik + 2 iff 



+ €S(v-n) = 0, (lb) 



and in case B they are written as 

Here, the normalized coordinates are £ = z/zq and 
i] = x/xq, where z and x are, respectively, the distances 
in the propagation and transverse directions, zq = Kxq, 
K is the propagation constant of the corresponding con- 
tinuous mode, /? is the effective wave vector mismatch, 
8(rf) is the delta- function, and g — c^zq is the normalized 
coupling constant between adjacent waveguides in the 
array, where ca is the coupling constant in physical 
units, as predicted by the coupled-mode theory. |T^ 
The normalized field amplitudes at the FF, <f> n and <E>, 
and at the SH, %p n and *, are proportional to their 
counterparts, u n ,U and v n ,V , measured in physical 




FIG. 1: Schematic of the proposed setting 
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FIG. 2: Field profiles of composite solitons. Here and in the 
next figure, top and bottom panels correspond to the solitons 
of types A and B, respectively, whereas left and right panels 
display odd and even solitons. Vertical lines designate the 
location of discrete waveguides. The wave vector mismatch is 
/3 — 0, and the soliton wavenumber is A = 3 and A = 1.5, in 
cases A and B, respectively. 



units: <pn = {zp/ ' y/xa)^2j<Tfc Un, \E = z j d V and ip n = 
(z /x )j c v n , $ = (z /^/xE)^2-fdlcU„U, where -y d ,c = 



led = 9P J A dAe* c (2w) •x(w,w)e n , c (w)e n , c (o;) 

in cases A and B, respectively, with P d and P c the 
power in the stripe waveguide and the power density in 
the slab waveguide, Ad and A c the transverse areas of 
the stripe waveguide and slab waveguide between two 
adjacent stripes, respectively, and e n (x,y) and e c (;r) 
fundamental modes of the stripe and slab waveguides. 
In this paper, we focus on the case of g — 1, which 
adequately represents the generic situation. 

The above equations were derived for the case when 
the FF and SH fields have the same (TE or TM) polar- 
ization. Using the general coupled- mode description |l8j , 
the system can be extended for the case of two polariza- 
tions, which will lead to a type-II x' 2 ' interaction^ Q 
involving two different FF components. This case will be 
presented elsewhere. 

Both variants of the model (A and B) neglect direct 
FF-SH interactions in the same waveguide, as they are 
suppressed by the large natural mismatch, while we as- 
sume that care is taken to minimize the mode mismatch, 
between the continuous and discrete waveguides. Such 
requirements can be easily fulfilled, as the geometry gives 
rise to different propagation constants for the same fre- 
quency in the waveguides of the two types. A model 
which takes into regard the residual SH-FF coupling in 
each waveguide can be easily considered too, but cases A 
and B, as defined above, are the most interesting ones. 
Note also that in this geometry the overlap between the 
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FIG. 3: The total power P of the semi-discrete solitons vs. 
the wavenumber A. 



FF and SH is smaller, as compared to the case in which 
FF-SH interactions in the same waveguide are employed. 

Composite solitons amount to stationary solutions of 
systems ([Ql and (J2J) in the form of c/>„(C) = </>n exp(iA£), 
tffo.C) =_*(?7)exp(2iAC) and ip n (C) = i> n exp(2i\(), 
®{V> C) — ®(v) ex P(*^C)j respectively, where A is the soli- 
ton's wavenumber. Inserting these expressions in the un- 
derlying equations, we solved the resulting systems by 
dint of the Newton-Raphson method. Similar to the or- 
dinary discrete solitons, the composite ones can be odd or 
even: the former ones are centered at a site of a discrete 
waveguide, whereas even solitons are intersite-centered. 
Typical examples of odd and even composite solitons are 
displayed for both cases, A and B, in Fig. [3 

To examine the stability of the composite solitons, we 
have first applied the Vakhitov-Kolokolov (VK) criterion 
|l9j | , which predicts that the necessary stability condition 
is dP/dX > 0, where P is the total power [which is the 
single dynamical invariant of Eqs. and (J5J], P(X) = 

2 f+™ I {*. $ i I'*? + £„ I {0ni <M I 2 ' in cases A and B > 
respectively. The result is displayed in Fig. |3 which 

shows that (as may be expected) the solitons exist above 
the band of linear waves of the discrete subsystem, i.e., 
for A > 2 in case A, and A > 1 + (3/2 in case B, and 
both odd and even solitons are VK-stable in most of their 
existence domain. The prediction is significant, as even 
solitons are always unstable in ordinary discrete systems. 
Note that the power of odd solitons is slightly smaller 
than that of even ones. 

A noteworthy feature is that A-type solitons exist up 
to P = 0, while their B-type counterparts have a cutoff 
in terms of P. This can be explained by noting that the 
limit of P — > corresponds to the limit of broad (quasi- 
continuous) small-amplitude solitons. Straightforward 
analysis of Eqs. (|la|l demonstrates that this limit exists 
indeed, reducing to ifa + <f> m + |c/>| 2 <?V (2 (4 — (3)) — 0, 
* = 4> 2 / (2(4 -/?)), with {0,*} = {0„e^ 2 <,*e- 4 <}. 
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FIG. 4: Top panels: field profiles of twisted composite solitons 
of types A (left) and B (right). In both cases, the wave vector 
mismatch and soliton's wavenumber are j3 = and A = 3.6. 
The bottom panel shows the total power P vs. A, the inset 
blowing up a vicinity of the cutoff point (stability border) on 
the curve for /3 — 0. 

On the other hand, the same limit, if applied to Eqs. (|2a.p. 
amounts to a linear asymptotic equation, + = 0, 
which gives rise to no soliton solutions. 

Full stability of the solitons was examined in direct 
simulations of Eqs. Q an d ©? which completely corrob- 
orate the predictions of the VK criterion. In particular, 
those composite solitons of the B-type which are VK- 



unstable as per Fig. |3| decay into linear waves. 

We have also studied the twisted solitons, built as out- 
of-phase bound states of two odd ones, see examples of A- 
and B-type solitons with a twisted FF component in Fig. 
0] Their stability was also examined by means of both the 
VK criterion (see Fig. and in direct simulations. The 
results demonstrate that the semi-discrete twisted soli- 
tons exist if their wavenumber A exceeds a cut-off value, 
which depends on the mismatch /3, and they are stable 
in almost the entire existence domain. By contrast, in 
continuous media twisted solitons are always unsta- 
ble 0- EI ■ I n a small region near the cut-off the twisted 
composite solitons are unstable too. The cut-off being 
well separated from the band of linear waves, in simula- 
tions the unstable twisted solitons do no decay into linear 
waves, but rather evolve into a fundamental odd soliton 
or split in two such solitons. 

To summarize, we have shown that stable composite 
solitons, fundamental and twisted ones, can be supported 
by a slab waveguide coupled to an array of waveguides, 
both made of a material. This new species of semi- 
discrete solitons features unusual properties, viz., sta- 
bility of intersite-centered fundamental solitons, and of 
states with the twisted continuous component. One can 
envisage that this system may also support walking po| 
semi-discrete solitons, and that similar solitons can be 
found in dual discrete-continuous systems with the Kerr 
nonlinearity, as well as in multidimensional x systems. 
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